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We show using bosonization methods that extended Hubbard-Heisenberg models on two types of 
two leg ladders (without flux and with flux vr per plaquette) have commensurate pair-density wave 
(PDW) phases. In the case of the conventional (flux-less) ladder the PDW arises when certain filling 
fractions for which commensurability conditions are met. For the flux tt ladder the PDW phase is 
generally present. The PDW phase is characterized by a finite spin gap and a superconducting order 
parameter with a finite (commensurate in this case) wave vector and power-law superconducting 
correlations. In this phase the uniform superconducting order parameter, the 2kp charge-density- 
wave (CDW) order parameter and the spin-density- wave Neel order parameter exhibit short range 
(exponentially decaying) correlations. We discuss in detail the case in which the bonding band of 
the ladder is half filled for which the PDW phase appears even at weak coupling. The PDW phase 
is shown to be dual to a uniform superconducting (SC) phase with quasi long range order. By 
making use of bosonization and the renormalization group we determine the phase diagram of the 
spin-gapped regime and study the quantum phase transition. The phase boundary between PDW 
and the uniform SC ordered phases is found to be in the Ising universality class. We generalize 
the analysis to the case of other commensurate fillings of the bonding band, where we find higher 
order commensurate PDW states for which we determine the form of the effective bosonized field 
theory and discuss the phase diagram. We compare our results with recent findings in the Kondo- 
Heisenberg chain. We show that the formation of PDW order in the ladder embodies the notion of 
intertwined orders. 
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I. INTRODUCTION 

A spectacular dynamical layer decoupling of the trans- 
port properties accompanied by different types of order 
parameters (charge, spin and superconducting) develop- 
ing together has been observed experimentally in the 
stripe-ordered (or nearly ordered) cuprate superconduc- 
tor La2_a;Ba3;Cu04 at zero external magnetic fielcPl, m 
underdoped La2-2:Srj;Cu04 at moderate magnetic fields, 
and in optimally-doped La2-a;Ba2;Cu04 at low magnetic 
fields:- A sequence of phase transitions are seen in these 
materials with the "normal" to charge-stripe ordered 
transition occurring first, followed by a spin-stripe or- 
der transition with a lower critical temperature. For in- 
stance, in La2_a;Baj.Cu04 near doping x = 1/8 a spec- 
tacular decoupling of the layers develops in transport 
measurements with the ratio of the c-axis pc to the ab- 
plane pat resistivities, becoming larger than 10^, begins 
to develop quite rapidly at temperatures right below the 
spin-ordering transition. At a critical temperature of the 
order of T^^ ~ 20K (depending on the precise doping) 
the copper oxide planes appear to become superconduct- 
ing while the c-axis transport remains resistive. The full 
three-dimensional resistive transition is seen only below 
lOK. A superconducting state with a Meissner effect 
and (presumably) d-wave superconductivity is seen below 
T^^ - AK. However, even though t he critical tempera- 
ture of the uniform d-wave superconducting state is much 
lower near a; = 1/8 than for other doping levels (where it 
is typically ^ AOK), the experiments show that the anti- 
nodal superconducting gap is essentially unsuppressed.!^ 



A strikingly similar transport anisotropy has been ob- 
served very recently in the temperature-pressure phase 
diagram of the heavy fermion superconductor CeRhlus W 
In this strongly correlated material the orders that de- 
velop are conventionally identified as a spin-density wave 
metallic state and a uniform d-wave superconductor. In 
the phase in which both orders coexist the ratio pd Pab 
becomes large (~ 10^) with pab eventually becoming 
unmeasurably small as is the superconductivity became 
two-dimensional (as in the case of La2-2;Baa;Cu04 near 
1/8 doping!). 

The most unusual aspect of these experiments is not 
just the existence of multiple coexisting orders, but the 
dynamical layer decoupling seen in transport. That is, 
the existence of a significant temperature range over 
which there is a form of two-dimensional superconduc- 
tivity in the planes but which are otherwise decoupled as 
if there was no Josephson effect between them. 

Berg et aPti^l showed that these seemingly contradic- 
tory results can be explained if one assumes that in this 
state charge, spin and superconducting orders are not 
competing with each other but rather that they are in- 
tertwined, with the superconducting state in the planes 
also being striped, i.e. it is a unidirectional pair density 
wave (PDW) with the property that the phase of the su- 
perconducting order parameter alternates in sign, as if 
the axes of the d-wave order parameter were to rotate by 
90° (with vanishing average value for the superconduct- 
ing order parameter) . The order parameter for the PDW 
state is 

" ^ ' Wpdw — Wpdw ^ ' 
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with the PDW ordering wave vector Qpdw pointing in 
the direction normal to the stripes, and (in the LTT lat- 
tice structure of La2_a;Baa;Cu04) rotates by 90° from 
plane to plane. Translation invariance of the underlying 
system further dictates that the ordering wave vectors 
for spin (QsDw)i charge (Qcdw) and superconducting 
(QpDw) order parameters obey the relation 2(5pdw — 

2QsDW — QcDW- 

In this paper we show that PDW type states do oc- 
cur in the phase diagram of strongly correlated systems 
on two-leg ladders. We focus on these systems since in 
this context the physics of strong correlations can be con- 
trolled and there are powerful bosonization and numer- 
ical methods to investigate their phase diagrams and to 
compute their correlators. For reasons that will be ex- 
plained below, states of this type generally involve strong 
correlation physics which is much harder to control in 
two dimensions. Results obtained in one-dimensional sys- 
tems, but with lots of caveats, can be used to develop at 
least qualitatively theories of two-dimensionally ordered 
states. In this work we will use bosonization methods, 
which are accurate at relatively weak coupling, to show 
how intertwined orders of this type develop in two-leg 
ladder. This work in many ways is an extension of the 
results of Ref . JT] that showed that the PDW state rep- 
resents the spin-gapped phase ( "Kondo-singlet" ) of the 
Kondo-Heisenberg chain. 

Non-uniform superconducting states were proposed by 
Fulde and FerrelP^l and by Larkin and Ovchinnikov^, 
and are conventionally called FFLO states. The super- 
conducting component of the PDW state discussed above 
has the form of a Larkin-Ovchinnikov state.'i^lThe Fulde- 
Ferrell stat^^ has spiral order. FFLO states have been 
proposed over the years for different types of supercon- 
ducting materials, most recently in the phase diagram 
of the heavy fermion superconductor CeCoIns at finite 
magnetic fields,'^ although the experimental evidence for 
them is weak (at best). There are also recent theoretical 
proposals for FFLO-type states in cold atomic fermionic 
systems with unbalanced populations,^ and in quantum 
wires of multi-valley semiconductors.'^ 

In the conventional theory of FFLO states one assumes 
a BCS-type system (a Fermi liquid) in which the spin 
up and down Fermi surfaces are split by the Zeeman in- 
teraction with an external magnetic field. As a result 
the ordering wave vector of the FFLO states is tuned 
by the Zeeman interaction and, hence by the magnitude 
of the magnetic field. A problem with this mechanism 
is that the usual nesting of the Fermi surface that leads 
to a BCS state zero-momentum pairing is generally lost 
if the Fermi surfaces are split, and the superconduct- 
ing states with finite-momentum pairing can only hap- 
pen for large enough attractive interactions, instead of 
being an infinitesimal instability as in the conventional 
BCS case.l^ FFLO phases driven by the Zeeman inter- 
action (and hence with a spin- imbalance and a broken 
SU{2) spin symmetry) in quasi-one-dimensional systems, 
including ladders, were discussed in Ref. PJ]- 



On the other hand the PDW state is the result of 
strong correlation physics, does not involve a Zeeman 
interaction, and does not occur at weak coupling. In the 
perspective of Ref. 8 , the PDW state is another man- 
ifestation of the concept of frustrated phase separation 
which is behind the development of inhomogeneous elec- 
tronic states in strongly correlated materials,'^^! whose re- 
sult are electronic liquid crystal phase^^ of which the 
PDW is a particularly interesting example. An under- 
standing of the physics of this state should shed light 
on the connection between electronic inhomogeneity and 
superconductivity. 

A BCS-mean-field-theory of a extended Hubbard 
model in two dimensions has been developed by Loder 
and coworker ^ ' ^ ^ I who found that, as part of a rich phase 
diagram, the PDW is the ground state for large enough 
interactions. Earlier mean field theory results in the t — J 
model by Yang et aW^ found that the PDW is energeti- 
cally very close to being the ground state. These authors 
found that (within their mean field theory) the ground 
state is a modulated d-wave superconductor, i.e. a state 
in which the uniform d-wave order parameter coexists 
with the PDW state, i.e. a state in which supercon- 
ductivity and stripe order coexist .'SS Several variational 
Monte Carlo calculations in the 2D t — J model have 
found that this is a very competitive state near doping 
X = 1/8 although not quite the best variational ground 
state.'^SnlZl. While these results are encouraging they suf- 
fer from the problem that either they use approximations 
that are not controlled in these strongly coupled regimes 
(as in the mean field studies), or that the search yields 
the best variational state within a restricted class which 
would be adequate at weak coupling but not in these 
regimes. Nevertheless these results indicate that from 
the point of view of local energetics the PDW state is 
very competitive and would most likely be the ground 
state for some extension of the Hamiltonians that were 
studied. A very recent and numerically intensive calcu- 
lation using infinite projected entangled pair states (an 
extension of the density matrix renormalization group) 
by Corboz et aP^have found stripe states coexisting with 
superconductivity in the 2D i — J model but has not yet 
investigated the presence (or absence) of a PDW state. 

In this paper we revisit the problem of the phase dia- 
gram of extended Hubbard-Heisenberg type models on 
two-leg ladders using bosonization methods. Here we 
show that the two-leg ladder has a phase that can be 
identified with a version of the PDW state, another phase 
in which uniform superconducting order is present and 
there is a continuous phase transition between these two 
phases. The two-leg ladder is an ideal model-system to 
study since it is well known from DMRG and bosoniza- 
tion results to have broad regimes of coupling constants 
and doping in which the ground state has a spin-gap 
and c?-wave superconducting correlations, as well as a 
strong tendency to stripe order.'^^HH Using a weak cou- 
pling band terminology, the PDW state is stabilized when 
one of the band of the ladder, say the bonding band, is 



3 



at special commensurate fillings. 

While there are extensive boso nizati on studies of the 
phase diagram of two-leg ladders the existence of 
phases with PDW order has not been previously investi- 
gated. Here we show that there is a connection between 
two leg ladder and the Kondo-Heisenberg (KH) chain, a 
model of a ID electron system interacting with a spin-1 /2 
Heisenberg antiferromagnetic chain through Kondo inter- 
action. In a series of papers, Zachar et. aLl^ShsZI showed 
the existence of a phase in which the correlations of the 
uniform superconducting order parameter decays expo- 
nentially fast while the correlations of a staggered super- 
conducting state (composite SC) has quasi long range 
order. Berg et. alP^ studied the KH problem using the 
Density Matrix Renormalization Group (DMRG) calcu- 
lations and showed that the spin-gapped phase of the KH 
chairPSlis a commensurate PDW state with = A* -.. 
Similarly, the PDW phases that we find in the Hubbard- 
Heisenberg model are also commensurate and have wave 
vector QpDW = In addition, and similarly to what 
happens in the KH chain, the PDW order parameter is a 
composite operator of a triplet uniform SC order param- 
eter of the anti-bonding band and the antiferromagnetic 
(Neel) order parameter of the bonding band. Separately, 
these two order parameters have short range order in the 
PDW state. In the PDW phase in the ladder system 
translation symmetry is broken spontaneously whereas 
in the KH chain it is broken explicitly by the spacing 
between the static spins. We discuss in detail the quan- 
tum phase transition between the PDW state and the 
phase with uniform SC order and it is found to be in the 
universality class of the quantum Ising chain. 

We also consider an extended Hubbard-Heisenberg 
model on a two-leg ladder with flux $ per plaquette. 
An important difference of this ladder system is that for 
general flux time reversal symmetry is broken explicitly 
except for the spacial case of flux ^ — ir per plaque- 
tte which is time reversal invariant. The presence of a 
non- zero flux changes the band structure by doubling 
the number of Fermi points at which the bonding band 
crosses the Fermi energy. We have not explored in full 
the complex phase diagram of this system except for the 
case of flux $ = tt. Here we found that for generic fillings 
of the bonding band the system obeys an Umklapp con- 
dition which leads to a spin gap state. We explored the 
phase diagram in this case and found that it generally 
supports two types of uniform and PDW superconduct- 
ing orders. In this system the PDW order parameters 
are bilinears in fermion operators and are not compos- 
ite operators and in the previous case. Here, as in the 
conventional ladder and in the KH chain, there is no co- 
existence between PDW and uniform orders: when one 
order develops (which in ID means power law correla- 
tions) the correlators of the other order parameters de- 
cay exponentially at long distances. In addition to the SC 
phase we also found four incommensurate CDW phases. 
The quantum critical behavior of this system is more 
complex, reflecting the larger diversity of phases that we 



encountered. In particular while the generic quantum 
phase transitions are also in the universality class of the 
quantum ising chain, for some special choices of param- 
eters the symmetry associated with the quantum critical 
behavior is enlarged to U{1) and it is now in the univer- 
sality class of a spinless Luttinger model. 

Two-leg ladders with flux $ per plaquette were stud- 
ied (both analytically and using numerical DMRG meth- 
ods) by Roux and coworker^^ in their work on diamag- 
netic effe cts in two-leg ladders, as well as by Carr and 
coworker^^HEIl -^y^q uged bosonization methods to study 
many aspects of the phase diagram. However in their 
work these authors did not consider the case of flux $ = tt 
per plaquette in which time reversal invariance plays a 
key role, and the problem of PDW phases, that (as we 
show occur here) generically at flux $ = tt, and which is 
the focus of the present work. 

This paper is organized as follows. In section [H] we in- 
troduce the model of the two-leg ladder and its effective 
field theory using bosonization methods. We can draw 
the phase diagram using the microscopic parameters of 
the ladder in the the weak coupling regime where their 
relation with the coupling constants of the bosonized the- 
ory is known explicitly. Although the form of the effec- 
tive field theory does not change with the strength of 
the coupling constants, in more strongly coupled regimes 
numerical methods must be used to established this re- 
lation. In section |III| we present the bosonized theory 
of a ladder whose bonding band is half-filled. Here we 
show that this effective low energy theory has a hidden 
self-duality. In section |IV| we use renormalization group 
methods to determine the phase diagram for the case of 
a half-filled bonding band, and show that, in addition of 
a Luttinger Liquid type phase, it also has two SC states, 
one with uniform SC order and the other with PDW order 
with wave vector Qpdw = We show that in this case 
there is a direct quantum phase transition between the 
phase with uniform SC order and the PDW phase which 
is in the universality class of the quantum Ising chain.. 
In section |V] we extend this analysis to regimes with a 
bonding band at other commensurate fillings. The result- 
ing PDW phase with wave vector Qpdw coexists (or is 
intertwined) with an also commensurate charge-density- 
wave (CDW) state in the bonding band with wave vec- 
tor QcDW = QpDw/2. Unlike the half-filled case, this 
state does not occur at weak coupling. In section [VT| we 
consider an extended Hubbard-Heisenberg model on a 
two-leg ladder with fiux $ per plaquette. Here we show 
that the commensurate PDW phase arises naturally in 
this frustrated band structure although through a differ- 
ent mechanism. The conclusions are presented in section 
|VII[ The RG equations for the general case are presented 
in Appendix |A] The solution of the effective field theory 
of the flux $ = TT model for special combinations of cou- 
pling constants and refermionization is given in Appendix 

m 
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II. MODEL OF THE TWO-LEG LADDER AND 
EFFECTIVE FIELD THEORY 

Consider a model of the two-leg ladder whose Hamil- 
tonian is H = Hq + Hint ■ The kinetic energy term is 

-i^^{4j,.C2,J,a+h.C.} (2.1) 

with t and t± being, respectively, the intra-leg and inter- 
leg hopping amplitudes, and i = 1,2 being the chain 
index and j being the lattice site index . The interac- 
tion terms of the ladder Hamiltonian have the form of a 
extended Hubbard-Heisenberg model, 

i,3 i,3 

+Vj_ ^ nijn2,j + ^ {nijn2,j+i + ni j+in2,j} 

3 3 

i,3 3 

(2.2) 

where [/ is the on-site Hubbard repulsion, V|| , Vj_ and 
Vd are the nearest neighbor and next-nearest neighbor 
"Coulomb" repulsions, and Jy, Jj_ and are the nearest 
and next-nearest neighbor exchange interactions. 

In the weak coupling regime, [/, V, J ^ t^t^_^ we pro- 
ceed by first diagonalizing the kinetic energy term Hq 
and finding its low-energy spectrum. This can be done 
by switching to the bonding and anti-bonding basis de- 
fined as (at each rung j of the ladder and for each spin 
polarization a) 



Cb,a = ;^(^2 ± Ci) 



(2.3) 



In the new basis the kinetic term reads as 

7/=a,6 j,a 



in which b and a stand for bonding and anti-bonding, 

and where t^j = t ± tj_ for rj = b,a respectively. 

In order to find the continuum limit representing the 
low-energy and long-wavelength behavior of the model, 
we linearize the energy dispersion of each band of the 
ladder around the respective Fermi wave vector 



er,{k) K Ef+ Vr,{k - kF,n) 



(2.5) 




where w,, 2t^ sin(fc^^) are the Fermi velocities, kp ,,j the 
Fermi wave vectors for each band, and Ep is the Fermi 



FIG. 1. Schematic picture of the bonding and anti-bonding 

bands. The bonding band b is kept at half-filling. The filling 
of the anti-bonding band a is general. Here kp.a and kF,b = f 
are the Fermi points for each band. 



energy. We now consider the regime of small fluctuations 
close to the Fermi points of each band: 

^c^,,,. ^ i?^,<.(:r)e"=^''^ + L^,„{x)e-"'^'>^ , (2.6) 

V 0, 

where R and L are right- and left-moving components 
of the electron field, x = ja is the position, and a is 
the lattice constant (the rung spacing). In this limit, 
the kinetic term of the Hamiltonian takes the standard 
continuum form 

= X] / dx{-ivn) {Rl^^da,Rn,a - L\ ,^d^L^^„} . 

(2-7) 

The most general continuum interacting Hamiltonian 
density (up to possible Umklapp processes that will be 
discussed below) compatible with the charge conservation 
and the global SU (2) spin symmetry with scaling dimen- 
sion two (marginal) operators has the following form 

'Hint = X {/clrj {Jrti + Jlt)) + 9c1ti JrtiJlt]} 

rj—a,b 

+ E {ffT^vi-^Rri + ^Lri) + 9sri Jri] ■ Jli]^ 

+ fc2 {JRaJRb + JLaJhb) (2.8) 
+ gc2{Jn.aJLb + JhaJRh) 

+ fs2{JRa • JRb + Jha ' J Lb) 

+ 9s2(JRa ■ Jbb + Jba ' J Rb) 

+ \t{^l ■ A„ + h.c.) H- A,(A|;A„ -t- h.C.) 

in which Jr/l.ti £^re the right and left moving components 
of the charge density = Jr^j) + <^L,r) for each band rj 

jR,ri = '^Ri,Ti^a,r,^ '^L,r) = '^Ll,n^a,ri (2-9) 
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Jr/l.ti are right and left moving components of spin den- 
sity Jri — Jr^t] + Jl,7] for each band 77, 

(2.10) 

where the components of the vector t? are the three Pauh 
spin matrices, and 



. (2-11) 

are singlet and triplet pairing operators respectively for 
each band rj. 

In the weak coupling limit, the relation between the 
coupling constants of the continuum theory of Eq.(2.8) 



and the parameters of the Hamiltonian of the microscopic 
lattice model of Eq.(2.2) can be found through this naive 
continuum limit procedure. This has been done before 
and can be found, for example, in Ref. [35]. We note 
that here by keeping only (naively) marginal operators 
we have neglected a host of irrelevant operators that are 
present in the lattice model that do not change the form 
of the low energy theory (although change the definition 
of the coupling constants by small amounts) . In the inter- 
mediate to strong coupling limit, either non-perturbative 
methods such as Bethe ansatz (applicable only if the sys- 
tem is integrable) or numerical density-matrix renormal- 
ization group (DMRG) calculations are required to make 
a quantitative connection between the the lattice model 
and the effective continuum field theory that we will use 
below. Nevertheless the form of Hint is general as seen 
here. 

Here we will be interested in the case of a half-filled 
bonding band. This situation will happen naturally as 
the total filling of the ladder is varied without break- 
ing any symmetries of the ladder, we should note that 
if one wanted to specify the filling of the bonding and 
anti-bonding bands separately, it would be necessary to 
set a chemical potential difference between the two legs 
which would make the ladder asymmetric. At any rate, 
for a half-filled bonding band, in addition to the interac- 
tions presented in 'Hint, we need to include the following 
(bonding band) Umklapp process 

nu,b = guALl-^Ll^RbiRbte'^''^'" + h-c.) (2.12) 

The value oi gu,b in the weak coupling regime is found to 
be given by 

-9u.b = liU+V^)-iV\\+Vd)-lj± + ^{J\\+Jd) (2.13) 

We will furthermore assume that for a substantial range 
of parameters of interest Jiu is marginally relevant for a 
half filled bonding band. We will get back to this point 
in the next section. 



III. ANALYSIS OF A TWO-LEG LADDER 
SYSTEM WITH A HALF-FILLED BONDING 
BAND 

We will now consider in detail the case when the 
bonding band is half filled and its Fermi wave vector is 
kpb = /'^- In this case there is a marginally relevant 
interaction representing the Umklapp process mentioned 
above. The main effect of this process is to open a charge 
gap Ac in the bonding band. Therefore for the ener- 
gies much smaller than the charge gap, one can assume 
that the charge degrees of freedom on bonding-band b 
are frozen-out and hence play no roll in the low energy 
limit of the remaining degrees of freedom. Moreover, due 
to the charge gap in the bonding band, all the interac- 
tions with net charge transfer between the bands, namely 
singlet and triplet SC in Eq. (2.8) processes, become ir- 



relevant. Therefore the only remaining charge degree of 
freedom, which is that of the anti-bonding band a, is de- 
coupled from the rest of the dynamics, and this being a 
one-dimensional system, the effective field theory of the 
charge sector of the anti-bonding band is described by 
the Luttinger Liquid (LL) theory. In its bosonized form 
the effective Hamiltonian density for the charge sector 
involves the Bose field (j)c and its dual field 6c for the 
anti-bonding band a only (where to simplify the nota- 
tion we dropped the label) and reads reads as 



1 



{dx<f>c 



(3.1) 



where is the charge Luttinger parameter of the system 
and Vc is the velocity. The fields 0c and 9c are dual to 
each other and satisfy the standard equal-time canonical 
commutation relations 

[Mx),d,ec{x')] = iSix-x') (3.2) 

which identifies He — d^Oc with the canonical momen- 
tum. 

All the possible remaining interactions between the de- 
grees of freedom of both bands are in the spin sector. The 
most general interacting Hamiltonian for the spin sector, 
which is symmetric under the exchange of the band index 
1] — a,b, is 

"Hs = ~gsl{JRb ■ Jhb + -hia ' Jha) (3.3) 
— gs2{JRb- JLa + Jhb- Jro) (3.4) 

Following the standard bosonization procedure in one di- 
mension, in terms of spin boson fields 



± <t>s,a) 



we arrive at the following form for T-L^ 

cos(^/47^0s+) 
2(™)2 



(3.5) 



(3.6) 



.gsi cos(v47r0s_) +3^2 cos(v47r6'5_) 
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where the Luttinger parameters Ks± and velocities Vs± 
are related to gs± = (5^,1 ± gs.-i)l^ as 



K 



s± 



1 2TTVf + gs± 
Sttw/ - gs± 



Vs± 



(3-7) 

in which Vf is the Fermi velocity of the noninteracting 
problem. The dual fields {9s± , 4's± ) obey similar com- 
mutation relations as the dual fields in the charge sec- 
tor. In general gsi is different for each band gib ^ 
gia- This introduces terms involving operators of the 
form dx4>s+dx4's- and dxOs+dxOs-- Although these are 
marginal operators we will neglect them for now since 
we will later argue that the results are not essentially 
affected by these terms in phases with a spin gap. In 
the absence of the spin gap, i.e. in the Luttinger Liquid 
phase, these operators change (among other things) the 
scaling dimensions of the observables.H2 

Upon inspecting the Hamiltonian of the spin sector 
Eq. (3.6), we see that it is invariant under the duality 



transformation 



s-,ffsl,5s2,^s-) 



(3.8) 

Thus, this duality symmetry guarantees the existence of 
a dual phase associated with the vanishing of the cou- 
pling constant gs2 = 0. We will see that in contrast to 
the PDW phase which is controlled by KH fixed point, 
in the dual phase the uniform SC is the most dominant 
instability. We will discuss the implications of this sym- 
metry on the phase diagram later on. 

In the gsi — limit the Hamiltonian of Eq.(3.6) 
turns out to be the same as the effective field the- 
ory description of the continuum limit of the one- 
dimensional Kondo-Heisenberg (KH) chairP^IMlZl with 
nearest-neighbor Kondo spins. The KH chain is a system 
of a IDEG (which usually is taken to be non-interacting 
but this is not necessary) and a one-dimensional array of 
spin-1/2 degrees of freedom, a ID quantum antiferromag- 
net with exchange interaction Jh- The spacing between 
the spin degrees of freedom defines the unit cell of the 
chain and in general it is not equal to the lattice spacing 
of the IDEG. The coupling between the spin chain and 
the IDEG is the Kondo exchange coupling Jk- 

That these two problems have almost the same low 
energy effective field theory of the same form is not a 
coincidence since there is a formal analogy between the 
two problems. In both cases we have a gapless IDEG, 
the free fermion band of the KH problem with the IDEG 
of the electrons in the anti-bonding a-band of the two 
leg ladder system, which in both cases is coupled to a 
Heisenberg spin-1/2 chain. It is known that the KH chain 
(regardless of the spacing between the Kondo spins) has 
a broad regime of its phase diagram in which there is a 



spin gap.*^ This phase has been identified with a com- 
mensurate PDW phase.'^ One difference between these 
two systems is that in the two-leg ladder the coupling 
is the tunneling matrix element t±^ whereas in the KH 



case it is the local Kondo exchange Jk- However since 
the bonding band of the ladder has a charge gap the 
only active coupling allowed at low energies is also the 
effective exchange coupling. Thus in the g^i = limit 
both problems are the same. We will see below that in 
the regime with 5^1 = the parameters of the Hamil- 
tonian flow under the renormalization group to a stable 
fixed point of the characterized by pair density wave cor- 
relations. Also, in the g^i — limit the system has an 
exponentially small gap (which can be determined form 
a mapping to the SU{2) Thirring modeP^ that is sta- 
ble against small perturbations of the form we discussed 
abov(ii2j^. 

However, there is an important qualitative difference 
between these two systems. While the Kondo-Heisenberg 
chain is translationally invariant only if the lattice spac- 
ing of the quantum antiferromagnet (the distance be- 
tween the Kondo spins) is the same as the lattice spacing 
of the IDEG, whereas the ladder is a translationally in- 
variant system in all cases. This will play an important 
role in our analysis. 



IV. PHASE DIAGRAM OF THE SYSTEM WITH 
A HALF-FILLED BONDING BAND 



We will now discuss in detail the phase diagram 
and phase transitions of a extended Hubbard-Heisenberg 
model on a ladder in which one band, the bonding b band 
is half filled, as shown schematically in Figj2] 



A. Weak Coupling RG analysis 

The total effective Hamiltonian isH ~ Hc + Hs, where 
the Ha milt onian density for the charge sector T-lc is given 
by Eq.(3.1) and the Hamiltonian density of the spin sec- 
tor Hs is given by Eq.(3.6). The total Hamiltonian has 
five coupling/parameters, Kc, Ks±, and g± (or equiv- 
alently gsi and gsi)- The Luttinger parameter of the 



charge sector will not renormalize since the charge sector 
decouples, while all the parameters in the spin sector are 
subject to renormalization. The one-loop RG equations 
for the couplings in the spin sector are 



dKs+ 




dl 


87r2 


dKs- 
dl 


1 2 

- 8^2 5.2 


dgsi 
dl 


(2 - Ks, 


dgs2 
dl 


(2 - Ks, 



2 

2 ffsl' 



Ks 



-)9s 



(4.1a) 
(4.1b) 
(4.1c) 
(4. Id) 
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FIG. 2. Schematic phase diagram when the bonding band of 
the ladder is half-filled, shown as a projection of the SU{2)- 
invariant RG flows onto the {gsi,gs2) plane. The solid black 
lines represent the separatrix between different phases. Due 
to the duality symmetry in the effective Hamiltonian for the 
spin sector, the phase diagram is symmetric around gs- ~ 0. 
The quadrant gsi,gs2 > flows into the Gaussian fixed point 
gsi ~ gs2 = and is a Luttinger Liquid (LL) phase. Region 
II is controlled by the PDW strong coupling fixed point at 
(gsi = 0,gs2 = —oo). Region I is controlled by decoupled 
spin-gapped fixed point of {gsi = 0,gs2 ~ —<x). There is a 
KT transition from the LL behavior across the half-line gs2 = 
and gs.i > to region II, and across the half-line gai = 
and gis,2 > to region I. The half-line gsi = gs2 < represents 
a quantum phase transition between the two strong-coupling 
phases and it is in the Ising universality class. 



B. Luttinger Liquid Phase 

We start with the case where gsi = and |5s2(0)| small. 
The analysis for this regime is similar to what is done 
in the Kondo -Heisen berg chain (with nearest neighbor 
Kondo spins) .1^^1261371 



In this regime the RG equations of Eq. (4.1) become 



dl 
dgs2 _ 
dl 



dK, 



dl 



1 

8^ 



1 



(4.2) 



s+ 



)9. 



s2- 



To guarantee the SU(2) invariance, these equations are 
subject to the initial conditions Kg_{0) = K~_^{0) ~ 
1 — gs2(0)/47r. The relation Ks- = is an invariant 
of this RG flow. Upon implementing this constraint, the 
set of RG equations (4.2 ) can be further simplified to the 



Kosterlitz RG equations 



dx 
~dl 



1 

8^ 



2 

2 i's2i 



dgs 



dl 



2xgs 



(4.3) 



satisfy g = ±47ra;. Thus, a point on the SU(2) invariant 
trajectory x{0) = Kg-iO) - 1 ~ -5s2(0)/47r < flows 
into the Gaussian (free field) fixed point gsi — gs2 — 0, 
Ks+ — Ks- = 1 where the system is a (spin 1/2) Lut- 
tinger Liquid which is a gapless and hence scale invariant 
system. In other words, the fixed point of {gsi — 0, gs2 — 
oo) is perturbatively unstable and flows to the Luttinger 
liquid fixed point along the SU{2) invariant RG trajec- 
tory. At this Luttinger liquid fixed point all three dif- 
ferent components of the triplet SC order parameter 
on the anti-bonding a-band and the Spin Density Wave 
(SDW) order on both bands decay as a power law. This 
is the Luttinger Liquid phase, shown in Fig|2j and its 
correlators are standard. Since it has one gapless charge 
mode and two gapless spin modes it is a C1S2 Luttinger 
state (in the terminology of Balents and Fishei'^. 



C. The PDW phase 

In contrast, points along the other SU (2) invariant tra- 
jectory, x{0) = Ks-{Q) - 1 ~ -.g2(0)/47r > 0, flow to the 
strong coupling fixed point {gsi = 0,gs2 = — oo). At this 
fixed point both cos{VATr9s-) and cos(-\/47r^s+) acquire 
non- vanishing expectation values, the fields are pinned at 

values {6s-,4's+) = ^{Ns- , Ns+){wheTe Ns^± are both 
odd or even integers at the same time) and their quan- 
tum fluctuations are gapped. It is easy to see that the 
identity {cos y^9s-) = {cos y^(j)s+) holds along the 
SU{2) invariant trajectories. This observation will be 
useful later. 

This means that a extended-Hubbard Heisenberg 
model on a ladder when one of its bands is half-filled, 
has is a Kosterlitz- Thouless (KT) phase transition from 
the gapless Luttinger Liquid phase to the spin-gap phase 
described by the same strong coupling fixed point of 
the spin-gap phase ("Kondo singlet") of the Kondo- 
Heisenberg chain (with nearest-neighbor Kondo "impu- 
rities"). 

However, at this strong coupling fixed point observ- 
ables that involve the dual fields of either (j)s- and 9s+ 
have vanishing expectation values and their correlations 
decay exponentially fast. This results in short range cor- 
relations for the singlet uniform superconducting order 
parameter on the anti-bonding a band, A^, and of the 
Charge Density Wave (CDW) order parameter, Ocdw 
for both bands a and &, which have the charge-spin fac- 
torized form 



1 



A„, = ^e-'^2^«-» C0s{V2TT^s,a), 
1 



OcDW,a/6 = — e-'^2.^-/^ cos(v^0,,,/f,) (4.4) 



vra 



where x — Ks- — 1 <C 1. The solution of this flow equa- 
tion is well known. Here we will only be interested in 
the SU{2) invariant trajectories of the RG flow which 



where here a is the lattice spacing. For instance, in the 
case of the singlet uniform SC order parameter A^ for 
the anti-bonding a band, in the s± basis, its spin part is 
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decomposed as 

cos{V2Tr(j)s,a) =cos(\/7r(/)5_)cos(v/7r(/)3+) 

+ Sm{y/TT(l)s-)sin{y/TT4>s+) (4.5) 

The presence of cos(-y7r0s_) and sm{y/TT(j)s-) guarantees 
exponentially decaying correlation functions. A similar 
form is also found for the CDW operators of the bonding 
and anti-bonding bands. 

However this does not imply that there are no long 
range correlations at this fixed point since there are com- 
posite order parameters built from products of these op- 
erators which do have quasi long range ordei^SI. To this 
end we define the order parameters for the PDW phase as 
the staggered part of the following product of operators 
in the two bandJ^l 



= Aa-Sl, = Aa-Jb + i-ir/^'O 



PDW 



(4.6) 



where St ~ Jb + {—l)^^'^Nb, and Jb is the total spin 
density vector on the bonding b band and Nb is the Neel 
order parameter also for the bonding b band. The explicit 
bosonized expression for the pair-density-wave order pa- 
rameter OpDW is 



O 



PDW 



^Aa-Nb = 



1 



■ cos(v27r(/)c,6)e 



2(7ra)2 

2 cos(\/47r^?s_) + cos(\/47r(/)s_) — cos{V4Tr(f>s+) 

(4.7) 

It is easy to see that the spin part of Opow has a nonzero 
expectation value in this phase. Therefore in spit of the 
fact that the correlation functions of the individual Neel 
(or SDW) order parameter of the bonding b band and 
uniform triplet SC of the anti-bonding a band are expo- 
nentially decaying, the correlation function of their prod- 
uct, the PDW order parameter OpoWi exhibits power law 
correlations: 

(Opdw(2;)O;dw(0)> ^C.^C^Ixp^/^- (4.8) 



in which Cc — (cos(v27r(/)c,6)) and Cs = (cos(v47r6's_)). 
Therefore Apow has quasi long range order with the 
exponent of 1K~\. The very same argument holds for 
a composite operator obtained from the product of the 
CDW order parameter on the bonding b band and the 
uniform singlet SC on the anti-bonding a band. Similar 
to the OpDW case, in gs\ = limit the individual CDW 
on the bonding b band and uniform singlet SC on the 
anti-bonding a band both decay exponentially fast but 
their product has quasi long range order with the same 
exponent as PDW order parameter. This is the Region 
n phase shown in Figj2] 



D. Uniform Superconducting Phase 

Let us now consider the opposite limit of = 0. Here 
one can repeat a similar RG analysis as in the gs\ = 



limit (or make use of the duality symmetry hidden in the 
problem, shown in Eq.(3.8)) to show that the RG with 
< 5si(0) <C 1 will be rcnormalized back to g^i = 
along the SU{2) invariant trajectory and we are again 
in the Luttinger Liquid phase we found before. There- 
fore, just like the fixed point at (gsi — 0,gs2 — oo) in 
the former regime, the fixed point (gsi = oo,gs2 = 0) is 
also not accessible in this case. This means that compo- 
nents of the uniform triplet SC or the Neel (SDW) order 
parameters have power law correlations. 

However, assuming again SU{2) invariance, in a sys- 
tem with gs2 = a small negative gsi < will flow to 
the strong coupling fixed point {gsi = —00,3^2 = 0). 
This is the fixed point dual to the {gsi = 0,gs2 = —00) 
fixed point under the duality transformation of Eq. ( |3.8[ ) . 
In this phase now the expectation values {cos{y^4^4's±)) 
are nonzero and observables that are functions of (j)s± 
have long-ranged correlations. Instead the expectation 
values (cos(\/47r0s_)) = and its fluctuation has short- 
ranged correlations. At this strong coupling fixed point 
the semi-classical expectation values of the (j)s^a and 0s, 6 
are such that (cos(\/47r0s+)) x (cos(A/47r0s-)) > and 
hence ^/4Tr(j)s+ — \/47r0s- = 0, tt mod 2n. In a phase 
controlled by this fixed point the two sectors, ±, have 
separate spin gaps. Furthermore, in this regime the ex- 
pectation value (cos(-\/27r0s.a)) 7^ is nonzero, and there- 
fore in the phase controlled by this fixed point the uni- 
form singlet SC on the anti-bonding a band has quasi 
long range order. 



{A,X^)Al{0))^ClJx\~^/^^ 



(4.9) 



where Csa = (cos(v27r(/)s_a)) in this phase. By the same 
argument, the CDW order parameter of the anti-bonding 
a band has quasi long range order as well. 



(OSdw(^)O 



CDW (0)) 



c 



(4.10) 



However, given the repulsive nature of the interactions 
the Luttinger parameters obey K^^a < 1 and therefore SC 
is the dominant fiuctuations of this phase. Nevertheless, 
unlike the PDW phase, this phase has dominant uniform 
SC correlations. Nevertheless we note that in this phase 
there exists subdominant CDW correlations. 

On the other hand we will now see that the correla- 
tion function of OpDW decays exponentially fast in this 
phase. Similar to the previous discussion, it is easy to see 
that (cos(A/47r0s+)) = (cos(A/47r0s-)) for a SU(2) invari- 
ant RG flow. Therefore, looking back at the structure of 
the PDW order parameter given in Eq.(4.7), we see that 
the expectation values of cos(-\/47r0±) cancel each other 
out at this fixed point. Hence, the expectation value of 
the spin part of OpDW is zero since the expectation value 
{cos{^/4^T9s-)) — vanishes exactly in this phase. More- 
over the two-point correlation function of Opow has to be 
proportional to the vertex operator cos(-\/47r6's-) whose 
correlations decay exponentially fast. Consequently the 
correlations of OpDW are short-ranged. 
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On the other, since there are independent spin gaps on 
both bands, the product of SC on a-band and CDW on 
6-band still has long range order . Therefore this product 
has similar correlations at both phases. This means it can 
not be used as an order parameter to distinguish these 
two phases. Therefore OpDW is the unique order param- 
eter to distinguish the state with PDW order from the 
state with uniform SC order at the spin-gapped regime 
of the two-leg ladder system with one band kept at half- 
filling. This is the Region I phase shown in Fig j2] and it is 
identified with a phase with uniform superconductivity. 



E. The PDW-uniform SC Quantum Phase 
Transition 

We will now discuss the quantum phase transition be- 
tween the state with uniform SC order (with power law 
correlations) and the PDW state. To this end we first 
note that both PDW and uniform SC in their associate 
phases have quasi long range order with the same expo- 
nent of This happens since the exponents are 
controlled in both cases by the decoupled charge degree 
of freedom left in the system. However the two phases 
are distinguished by the fact that one SC state is stag- 
gered (the PDW state) while the other is uniform. There- 
fore by symmetry the phase transition between these two 
states is similar to the transition from a state with trans- 
lational symmetry to the one with only a Z2 discrete 
broken translational symmetry. So we expect to find that 
this transition to be in the Ising universality class. How- 
ever, as we will see, the way this happens is actually 
rather subtle. 

To discuss the nature of the phase transition between 
Region I and Region II of Figj2] we need to look at the 
<?si = gs2 line in the parameter space. The duality sym- 
metry of the Hamiltonian implies that the phase diagram 
must be symmetric under it. We will now see that there 
is a direct quantum phase transition at the self-dual line, 
the half-line that separates Region I from Region II in 
Figj2] From the first RG equation in Eq.(4.1 1, it is clear 
that the Luttinger parameter K^^ — > flows to zero 
whenever either gsi or is relevant. Furthermore along 
the Qsi — gs2 line system flows to a new strong coupling 



fixed point with (cos(v47r(/)s+)) ^ and the field (j)s+ is 
pinned. At the strong coupling fixed point the effective 
bosonized Hamiltonian can be further simplified to the 
following 



+ ^cos(V4^0,_) 



Me 



cos(V4^6i,_), (4.11) 



where the new couplings ^x^/g are related to the gs's as 

5s(l,2) 



2tt 



(cos(V47r03+)). 



(4.12) 



of Eq.(4.11| can be re-fermionized in terms of two sets 
of chiral Majorana fermions. For = /i^, i.e. along 
the self-dual line of our problem, one of the chiral Majo- 
rana fermion pairs becomes massless and its mass changes 
sign accross this phase transition. Therefore, the phase 
boundary between the phases corresponding to gsi and 
gs2 is the in the same universality class as the quantum 
Ising chain, a theory of a (non-chiral) massless Majo- 
rana fermion. We should point out that the expression 
of the Ising order (and disorder) operators in terms of 
the bosonized fields of the ladder is highly non-local. 

This question can also be addressed directly through 
RG equations as well. Indeed along the self-dual line they 
reduce to 



2 7^2 2 



dl ^ K, 



(4.13) 



This system was discussed extensively by Lecheminant 
et. alP^, who showed at Ks- = 1 the resulting Hos 



Starting with Ks-{0) = 1 and \n0iO)\ > |m0(O)|, we 
see that RG flows to the the fixed point of (|/ie| — 
00, /i0 — 0,Ks^ — 00) which is the same as the PDW 
fixed point. In contrast, if we start with Ks-{0) = 1 
and 1/^^(0) I > |/ie(0)|, the RG flow will take us to 
{pe — 0, l/i^l — 00^ Kg^ = 0), i.e. the uniform SC fixed 
point. The flow with the initial condition ifs- (0) = 1 and 
9si{0) — 3s2(0) < will flow to the Ising critical point 
with gsi = gs2 = -00 and Ks- = 1 while 5^1 = 3^2 > 
flows to the Gaussian fixed point. 

Figj2] shows the results of a numerical calculation of 
all the S'?7(2)-invariant RG flows projected onto the 
(<?si,5s2) plane. Solid black lines represent the separa- 
trix between different phases. The low energy behav- 
ior of all the models with 5J7(2)-invariance which sat- 
isfy gsi,9s2 > is controlled by the Gaussian fixed 
point. The rest of the flow, except for the semi line 
.gsi = gs2 < 0, wiU end either at {gsi,gs2) = (-00, 0) or 
a-t (5s1j5s2) — (0, —00) depending on the initial values of 
the couplings. The RG flow pattern is symmetric around 
the gsi — gs2 line as dictated by the duality symmetry. 



V. BONDING BAND AT GENERAL 
COMMENSURATE FILLINGS 

We will now discuss the case of a ladder with a bonding 
band at other commensurate fillings. To understand this 
case it is useful to recall first the physics of the simpler 
problem of the extended repulsive one-band Hubbard- 
Model at quarter filling kp — n/A. This system has a 
quantum phase transition at a critical value of the near- 
est neighbor "Coulomb" interaction between a Luttinger 
liquid and a commensurate insulating CDW state. An 
intuitive classical picture of the ground state of such a 
system is as if the electrons occupy every other site on the 
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lattice with their spins arranged in a "stretched" Neel an- 
tiferromagnetic state, as in an antiferromagnetic Heisen- 
berg chain with twice the lattice constant. In this regime, 
the charge sector of the Hubbard chain is gapped (and 
hence insulating). 

This phase transition is driven by a higher order Umk- 
lapp interaction, which appears at third order in pertur- 
bation theory, which stabilizes this period 2 CDW state. 
As it is well known, and easy to see, in the bosonized for 
of this system the Umklapp term for the 1/4 filled band 
has the form (see, e.g. Ref. 3?] and references therein) 

■H„,i/4 = .91/4 cos(4%/27r(/)c) (5.1) 

The scaling dimension of this Umklapp operator is 'iK^^ 
where is the charge Luttinger parameter of the ex- 
tended Hubbard chain. Therefore, this Umklapp process 
is relevant for Kc < 1/4 which always lays in the inter- 
mediate to strong coupling repulsive regime. Although 
in this regime the bosonization formulas that relate the 
parameters of the microscopic model (in its naive con- 
tinuum limit) and the bosonized theory is no longer ac- 
curate, the form of the effective low energy bosonized 
theory retains its form as it is dictated by symmetry. 
The main problem is that the connections between the 
Luttinger parameter(s) and the microscopic parameters 
is more complex due to finite renormalizations induced 
by the irrelevant operators. In practice this relation must 
(and is) determined from numerical calculations on the 
microscopic model. 

For a system on a two-leg ladder one can pursue a 
similar line of reasoning and use the fact that, for cer- 
tain fillings, there is a similar Umklapp processes for 
the bonding band (for instance). However, just as in 
the case of the extended Hubbard chain, here too the 
couplings corresponding to these Umklapp terms in the 
strong coupling effective theory can not be easily re- 
lated to the microscopic parameters of the original lattice 
model and requires serious numerical work. Nevertheless, 
such Umklapp processes still exist and should eventu- 
ally become relevant. At this value of the parameters, 
where = 1/4, the Umklapp process for the bond- 
ing band becomes marginally relevant, and the system 
has a Kosterlitz-Thouless phase transition to a period 2 
commensurate CDW state that coexists with antiferro- 
magnetic order, with power law correlations in the spin 
sector. 

We can now follow the same approach we used for 
the case of a half-filled bonding band of the preceding 
sections to determine the phase diagram. We will not 
present the details here but give the main results. As 
before, the phase diagram in general has three phases: 
a) a Luttinger Liquid phase (similar to the one discussed 
by Wu et aP'^), b) a phase with uniform superconduct- 
ing order (and hence a spin gap), and c) a PDW phase 
(also with a spin gap). However, the ordering wave vec- 
tor of the PDW state is now Qpdw — 7''/2 (instead of 
QpDW — for the case discussed in the preceding sec- 
tions). This PDW state there is a "composite" CDW 



quasi ordered state (with power-law correlations) with 
degrees of freedom partially on the anti-bonding band, 
in this case with wave vector Qc^-w = 'i'^p.a + Opdw- 
Thus the resulting PDW phase is similar to the one dis- 
cussed in Refs. |lli and i37]. In contrast, the phase with 
uniform SC order has a CDW state that develops on the 
anti-bonding bans alone and has a conventional 2kp a or- 
dering wave vector. In spite of these differences with the 
case of the half-filled bonding band, the quantum phase 
transition between the PDW phase and the phase with 
uniform SC order for the quarter filled bonding band is 
also in the Ising universality class. 

A state with very similar properties was found in the 
Kondo-Heisenberg chain for a Kondo lattice with pe- 
riod 2 (see Ref.|ll|). However, while in the KH chain 
translation invariance is broken explicitly by the assumed 
fixed spacing of the Kondo spins, in the case of the 
ladder translation invariance is broken spontaneously at 
the Kosterlitz-Thouless quantum phase transition we just 
discussed. However the spin gap (and the PDW state as 
well) can only develop once this CDW state is formed. In 
this sense this is an example of intertwined (as opposed 
to competing) orders in the sense discussed by Berg et 
alW 

This line of argument can, naturally, also be extended 
to states in which the bonding band has other fillings, 
such as 1/2" for n = 1, 2, and consider Umklapp terms 
of the form cos(2n-\/27r0c)- Although such terms will 
generally be present, their effects become relevant only 
for < 1 /n^ which lays deep in in the extreme repulsive 
regime for n > 2. Thus, unless the system has substan- 
tial interactions beyond nearest neighbors, the resulting 
higher order commensurate CDW and PDW states will 
be quite weak and difficult to stabilize. 



VI. PDW STATE IN AN EXTENDED 
HUBBARD-HEISENBERG MODEL ON A 
TWO-LEG LADDER WITH $ FLUX PER 
PLAQUETTE 

We will now introduce and investigate another lad- 
der model in which we can show has a PDW phase. 
More specifically we will consider an extended Hubbard- 
Heisenberg model in a two-leg ladder with flux $ per pla- 
quette (in units of the flux quantum (jjQ = hc/e). As usual 
the flux is introduced by a Peierls substitution which here 
reduces to assigning a phase ±$/2 to the electron hop- 
ping matrix elements along the two legs which now be- 
come complex and are complex conjugate of each other, 
ti — t2, where ii^2 are the hopping amplitudes on top 
and bottom leg (see Figjs] a) . In addition to the the hop- 
ping along the rungs, we assume a real hoping amplitude 
between the legs. 
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The free part of the Hamihonian of this system is 

Ho--tJ2 {^'^^'cl+lC^, + e-'^^/'cl+.c,^, + h.c 
j 

-*^E(^Ij^2,, + 4,,ci.,) (6.1) 



in which « = 1, 2 is the chain index referring to the top 
and bottom chains respectively and j € Z denotes the 
lattice sites. To the best of our knowledge this electroni- 
cally frustrated system has not been discussed previously. 
The interaction terms, iJint that we will consider are the 



same as in the conventional ladder system and are given 
in Eq.(E2l). 



We will see here that this model has a very rich phase 
diagram which we will not explore completely. However 
wc will show that PDW phases occur naturally although 
through a rather different mechanism that we found in 
the conventional ladder discussed in the previous sec- 
tions. We will discuss this problem using the same meth- 
ods bosonization we used above. 

We begin by constructing an effective field theory. In 
momentum space the free fermion part of the Hamilto- 
nian becomes 



" dk 



- 2tcos(fc + $/2)cl(fc)ci(/c) - 2icos(fc - $/2)4(fc)c2(A:) - tj_{cl{k)c2{k) + cl{k)ci{k))\ (6.2) 




FIG. 3. a) A two-leg ladder with flux "I> (in units of the flux 
quantum) per plaquette. Here 1 and 2 label the two legs, j and 
j + 1 label two consecutive rungs. The hopping amplitudes 
are t\ = t e'*^^, t2 = t e"'*''^ and t±. b) Schematic plot 
of the dispersion relation of the bonding band Eb{k) for a 
non-vanishing flux per plaquette, $ 7^ 0. 

For t± = the band structure consists of right and left 
branches centered at A: = $/2 and k = — $/2 respectively. 
For tj_ j^Q a, full gap opens up at the crossing point of the 
two bands and bonding and anti-bonding bands form. 

The Hamiltonian is diagonalized if we switch to the 
new basis defined by the orthogonal transformations 
Crj{k) — J2i ^'^rj i{k)ci{k) with the transformation matrix 
M{k): 



Mr,4k) = 



/cos(e(fc)/2) 
Uin(^(fc)/2) 



-sin(e(fc)/2)\ 
cos(C(fc)/2) ) 



(6.3) 



where rj = a,b stands for the bonding and anti-bonding 
bands, i — 1,2 labels the legs of the ladder, and ^ is 
defined as 



sin(^(fc)/2) = 



cos(e(fc)/2) = 



L{k) 



v/rT^?(fc) 

1 

y/l+U^k) 



(6.4) 



in which 



tj_u{k) = 2tsin($/2)sin(fc)+Y/(2tsin($/2) sm{k))'^ + t']_. 
The inverse transformations are defined by the inverse 



matrix ^ as following 



dk 

-.^^ 
dk 

2^' 



cos(^)c^(fc)+sin(^)c^(fc) 



-sm(^)c^(fc)+cos(^)c^(fc) 



2 



(6.6) 



where b and a label the bonding and the anti-bonding 
bands respectively. 

The dispersion relations for the bonding and anti- 
bonding bands are 



E,,{k) = -2tcos($/2)cos(fc)±y^(2tsin($/2)sin(fc))2 +t2^. 

(6.7) 

Band structures of this type appear in quantum wires 
with two pockets^^ and in ID electronic systems with 
spin-orbit interactions (with the leg label playing the 
role of the electron spin) (see, e.g. Ref . [47] and refer- 
ences therein). The dispersion relations Ejf{k) satisfies 
the symmetries 

E^{-k) ^ E^{k), Ea{<P + 27:, k)^-Eb{<P,k) (6.8) 



12 



For a wide range of parameters, t±/t and flux the 
bonding band lias the form sketched in FigjS] with two 
minima in momentum space. For the rest of this section 
we will focus on the regime of the parameters in which 
the Fermi energy lies below the hybridization gap of the 
bonding band, Ep < -Bh(O) -2icos($/2) - tj_. In 
this regime, the Fermi energy crosses the bonding band 
at four distinct points, ±fci,±fc2, while the anti-bonding 
band is empty, as shown in Fig. [3]b. 

Here we will consider an extended Hubbard-Heisenberg 
model on a ladder with flux <i> per plaquette. We will see 
now that this system has an interesting phase structure. 
We will analyze this system using the same bosonization 
methods as in the conventional ladder. For reasons that 
we will explain below we will focus on the spacial case of 
flux $ = TT per plaquette. 



A. Low energy continuum limit 

In this work we are interested in the regime in which 
the Fermi energy crosses only the bonding band and 
hence the anti-bonding band is empty. Furthermore we 
assume that the interactions are weak enough that we can 
focus only in the low energy excitations near the four 
Fermi points ±fci,±fc2 where the Fermi energy crosses 
the bonding band. Furthermore, for the more interesting 
case of flux $ = tt, the Fermi points obey the commen- 
surability condition fci + ^2 = vr. In this case we also 
have u(fci) — u{k2). Hence the parameter ^(fc) obeys the 
same identity and it is the same for both Fermi points. 
Henceforth it will be denoted by ^. 

By looking only at the low energy fluctuations around 
the Fermi points in the bonding band, the expansion of 
Eq.(6.6) reduces to the operator identifications 



1 



-.c^Aj) ^ sin(^)Li,(a;)e''=i- +cos(|)i?i,(x)e-"=^^ 



'2' 



+ sin(|)i?2a(a;)e"=^^ -f cos(|)L2.e-"=^- 

(6.9a) 

^C2,.(j) ^ cos(|)Li,(a;)e''=i^ + sin(|)i?i,(x)e-"=^^ 

+ cos(|)i?2.(x)e'^-=- + sin(|)L2.(a;)e-"=^- 

(6.9b) 

where we have used cos(^(— fc)/2) = sin(^(A:)/2) and 
^(/ci) — ^(fe) = which are both true for $ = tt, and 
where we have also projected out the anti- bonding band. 
We will treat the Fermi point labels as a flavor index 
/- 1,2. 

By inspection of the free fermion lattice Hamiltonian 
one can see that the Fermi momenta ki and k2 are essen- 
tially determined by the flux $ and by the filling fraction 
of the bonding band. In what follows we will ignore the 
contribution of the anti-bonding band since it is empty 



and its excitations have energies larger than the cutoff of 
the effective field theory. 

Following a similar discussion as in section [Hj the non- 
interacting continuum Hamiltonian becomes 



/ = 1,2 



where vi = — -gpUi and V2 = ~dk\t^2 ^-''^ the Fermi veloc- 
ities associated with the two Fermi points. For general 
flux $ there is no symmetry relation the Fermi points 
and the two Fermi velocities are different, vi 7^ U2. 

However, the the case of flux $ = tt per plaquette the 
energy bands have the additional symmetry fc — tt — fc. 
This symmetry reflects that fact that an exchange of the 
two legs, 1 -;-> 2, is in general equivalent to the reversal 
of the flux $ o — <& which is the time-reversed state. 
However due to the flux quantization, the states with 
$ = TT and $ = — TT are equivalent since the Hamiltonian 
is a periodic function of the flux with period 27r (corre- 
sponding to unobservable changes by an integer multiple 
of the flux quantum). On the other hand, from Eq.(6.7l, 



we see that for flux $ = tt the dispersion relations are 
also invariant under fc — > tt — fc which amounts to an ex- 
change of the two fermi points. Thus, in the case of flux 
$ = TT which insures that the Fermi velocities are equal, 
Vi = V2, for all fillings of the bonding band (and of the 
anti-bonding band as well). Therefore, for flux $ = tt, 
the symmetry of exchanging the two legs implies that the 
effective low energy theory must have a symmetry under 
the exchange of the flavor labels 1 and 2 (together with 
a chiral transformation which exchanges right and left 
movers) . 

In order to introduce all possible four-fermion intra- 
band and inter-band interactions, one considers an ex- 
tended Hubbard-Heisenberg type lattice problem, just as 
what we did for the two-leg ladder system of Section [Hj 
and construct the continuum theory for the present case. 
All four-fcrmion interactions for this system can be rep- 
resented by simple diagrams similar of the type shown 
in Fig|4] All the interactions can again be classified into 
charge and spin current- current interactions, singlet and 
triplet SC couplings or Umklapp processes with different 
commensurabilities . 

This means that the effective field theory for the 
present system has the same field theoretical form as 
the Hamiltonian of the two-leg ladder system given by 
Eq. (2.8 1. The only difference is that in the present 
case, the two sets of right- and left-moving labeled by 
the flavor index / = 1, 2 are low energy fluctuations of 
the bonding band. Moreover the connection between the 
couplings in the effective theory and the microscopic pa- 
rameters of the original lattice problem is different for the 
two problems. The two top diagrams in Fig. [4] represent 
singlet and triplet SC interactions between the 1 and 2 
while the lower diagrams corresponds to the most rele- 
vant Q ~ 2(fci -I- fc2) Umklapp processes. We will further 
assume that the Fermi points ±fci, ±k2 are such that no 
other Umklapp processes are allowed. 
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FIG. 4. Schematic representation of all the processes leading 
to uniform SC couplings and Q — 2(fci + ^2) Umklapp pro- 
cesses. The sum of the top two diagrams represent the uni- 
form singlet and triplet SC interactions while the two lower 
diagrams correspond to the Umklapp process. 



The discussion of the phase diagram of this system 
in the incommensurate regime is analogous to what has 
been discussed in the conventional ladder by many au- 
thors. C. Wu et. al.l^find that the only SC state in the 
phase diagram of the system away from the half-filling 
is the uniform s- or d-wave SC. Similar conclusions hold 
for the incommensurate regime of the this model which is 
in Luttinger liquid phase with two gapless charge modes 
and two gapless spin modes, C2S2. The only difference 
is that in this case these modes originate entirely from 
the bonding band. 

For general flux $ and for certain filling fractions of 
the bonding band Umklapp process involving separately 
the pairs of Fermi points ±fci and ±^2 become allowed. 
The physics that follows in these cases is similar to what 
we discussed for the conventional ladder in SectionllTland 
will not be repeated here. 

However, for flux $ 7^ a new type of Umklapp pro- 
cess, shown in Fig|4] becomes allowed. For this process 
to be possible Fermi momenta must satisfy the condi- 
tion Q = 2(fci + k2). This Umklapp process leads to the 
following interactions: 



A„3 «i"2 + ^u4 n[ 



h.C. 



(6.11) 



where n/ (with / = 1, 2) are the 2kp CDW order param- 
eters associated with the Fermi points at ztki and ±fc2, 
with ordering wave vectors 2ki and 2^2 respectively. Sim- 
ilarly, n/ is the associated SDW order parameters with 
the same ordering wave vectors. When the commensura- 
bility condition is satisfied this process is marginal and 
needs to be included in the effective low energy theory. 



However the commensurability condition ki + k2 = t: 
can only be met if the flux is $ = tt. Furthermore, in 
this case the system is commensurate for all fillings of 
the bonding band. For $ = tt the one-particle spectrum 

is given by Et{k) = — ^ At'^ siv? (fc) +t\ which satisfies 
Eii^TT — k) = Eb{k). Therefore if fci is a Fermi momentum 
so is fc2 = Ti'—fci- Hence for flux $ = tt the system remains 
commensurate for all electron fillings. We will see below 
that for $ = TT the pair-density-wave state exists for all 
values of the filling (with the Fermi energy in the bonding 
band). From now on we will restrict ourselves to the case 

of flux (f> TT. 

The bosonized Hamiltonian for flux $ = tt is (including 
the Umklapp process) 

n ^^-f {K^dj^+r + i^-+i(9,0,+)2} 



g^i cos(v47r(/)5_) +.9^2 cos(V47r6's- 

35 cos(V47r6'c-) + gu5 cos(V47r0c-) 

53 cos(^/47^^^s-) + 34 cos{V4t:4>s-) 

gu3 cos(%/47r6's-) + gu4 cos(V47r0s_) 

(6.12) 



+ 



cos(v47r(/)s+) 
2(7ra)2 

cos{V4Tr(l)s+) 
2(7ra)2 

cos(V47r6'c-) 
2(7ra)2 

cos(-\/47r(/)c_) 
2(7ra)2 



where (j)± = {(f>2 ± 4'i)/V2 and similarly for 6 fields. As 
before, there are marginal operators (both in the charge 
and spin sectors) of the form dx4)+dx4>- and dxO+dxO-- 
However, as in Section [Hj these operators can be ignored 
since their main effect is a renormalization of the scaling 
dimension^^ which here translate in smooth changes of 
the phase diagrams (without changing their topology) 
and in the spin gap phases they have essentially no effect. 
The first two lines of Eq.(6.12) is the sum of four 



different LL Hamiltonians for for the two charge and 
spin sectors. The third line corresponds to different spin 
backscattering processes, while the fourth and fifth lines 
represent the singlet and triplet SC couplings and the 
Q = 2{ki + k2) = 27r Umklapp processes respectively. In 
addition to the relation between initial value of the lut- 
tinger parameters of different sectors and the couplings in 
the various current-current interaction given in Eq. (3.7 1, 
the spin SU(2) invariance dictates that 55 = 34 + 33 and 
9u5 = 5ti4 + 3ti3- This will be useful in the discussion of 
the RG equations and phase diagram. 



The Hamiltonian of Eq.(6.12| has several symmetries 



Similarly to the half-filled bonding band case discussed 
in the preceding sections, we find a duality symmetry 
in the s— spin sector, {(j)s-,Os-) — > {0s-, —4>s-), under 
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which the Hamiltonian of Eq.(6.12| retains its form. We 
will denote this symmetry by . However, self-duality 
holds only if .g^i = 5^2, 93 = 94 and gus = g„4- In ad- 



dition, the last two lines of the Hamiltonian of Eq.(6.12) 
have identical form which indicates that we can define 
yet another duality symmetry of the same form but this 
time in the c— charge sector, ((/)(,-, 6*0-) — ^ i(^c-,—4>c-)^ 
and which will be denoted by Zj". Self-duality in this 
sector requires, in a dditio n, to set 55 = 5,(5. Finally, the 



Hamiltonian of Eq.(6.12) is also even in the fields 4'c,±, 
9c,±, 4's,± and 0s,±, which reflects the invariance under 
the exchange of the labels of the Fermi points (or flavors) , 
1 O 2, which is an exact symmetry only for flux $ = tt. 

In the next section we will look at the different SC 
and CDW states, each with a unique symmetry proper- 
ties under the action of the total symmetry group, and 
construct the order parameters for each state in order to 
identify the associated quantum phase diagram. 



B. Order parameters and phases 

To identify all the phases present in the phase diagram 
we construct the associated order parameters consistent 
with the symmetries of the current problem. In terms of 
the two flavors / — 1,2 of the fermions in the bonding 
band we define for this system two uniform SC order 
parameters, A±, and two PDW order parameters A-t 
(both with ordering wave vector Qpdw — t^)- They are 

A± ~ {L2^Rii + Ri^L2i) ± {Li^R2i + i?2tii^) (6.13) 

Similarly we also four CDW order parameters, n± and 
n±, 



n± =^(iLi?i,±LLi?2. 



(6.14) 



and their adjoint operators. 

The relation between these order parameters and the 
microscopic pair fields and CDW fields is as follows. 

The pair fields defined on site j on each leg i — 1,2 
of the ladder. A*, on the rung j, Aj^, and on each leg 
i = 1,2, AJ are defined by 

A,;.,- =c 



A,^ =Ci J- ^C2j,; + C2,j,tCl,j.i 



a;. =c,,^(j)c,,^(j + 1) + c,,^(j + l)c,Jj) (6.15) 

These observables can be written in terms of the slowly 
varying chiral Dirac fermions Rf,a and i/,o- (for the 
two flavors / = 1,2) in the symmetrized and anti- 
symmetrized forms (with respect to the exchange of the 



labels 1 and 2 of the legs of the ladder) 

A] + A2 ^ sine A+ + (-1)"/"A+ (6.16a) 

A] - A] -cose(-l)'^/''A_ (6.16b) 

Af ^ A+ + sine {-ly'"' A+ (6.16c) 

A]_^+i -f Aj^j+i ^ 2 sine sin(ga/4) A_ 

_ (_i)^^2icos(9a/4) A_ (6.16d) 

A],^.+i - A^^.+i ^ -(-l)-/'^2icose cos(qa/4) A+ 

(6.16e) 

where q = 2{k2 — fci) and where we have used the defi- 
nitions of Eq.(6.13). We see that the SC order parame- 



ters A± and A± represent two different types of uniform 
SC states and PDW SC states (both with wave vector 
Qpdw = t^) respectively. These pairs of SC states dif- 
fer by their symmetry transformations under flavor ex- 
change. It is worth to note that in the flux <I> = tt 
model the PDW order parameters are actually bilinears 
of fermion operators, c./. Eq.(6.16|. This is in contrast 



to wha t w e found in the conventional two-leg ladder in 
section HI and to the recent results by Berg et aP^ in the 
Kondo-Heisenberg chain, where the PDW order param- 
eter is microscopically quartic in fermion operators. In 
this sense the PDW states of the flux $ = tt two-leg lad- 
der is closer in spi rit to the conventional construction of 
FFLO states,E^'^ven though the spin SU{2) symmetry 
is preserved here and explicitly broken in the standard 
FFLO construction. 

Similarly we can relate the site Ui j (with i = 1, 2 the 
leg index) and rung, nj^ electron charge density operators 



,12 



2lt 



(6.17) 

which, after symmetrizing and anti-symmetrizing with 
respect to the exchange of the two legs lead to a set of 
four order CDW parameters, n± and n±, 

The relation between the microscopic charge density 
operators of Eq.( 6.17 1 and the slowly varying chiral Dirac 



fermions i?/,cr and i/,cr (with / = 1,2) is 
ni^j + n2j 



>J-?+J2° + 



ni. 



"•2, 



nf - nf 



sine {-If'" e'«^/2n+ + e'«^/2^+ + h.c. 

(6.18a) 

- cose(ji - jl) - (cose e'''^/^ri_ + h.c. 

(6.18b) 

-^sine(j2+J-?) 

+ (-1)^/'^ e'«^/2„+ + sine e'«^/2^+ + h.c. 

(6.18c) 



-cose {-ly'" e''^/2 



h.c. 



(6.18d) 



where we used the definitions of Eq.(6.14), and the usual 



definitions of the (normal ordered) currents and densities 
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of the Dirac fermions (again with / = 1,2) 



and n± are 



Jf 



(6.19) 



We can also define CDW order parameters on the legs 
of the ladder. However we will not discuss them since 
it turns out that they can also be expressed in terms of 
the same four slowly varying observables n± and n± and 
hence do not bring new information. 

From these results we see that in general we find both 
uniform SC order parameters and PDW order parame- 
ters, which always have a commensurate ordering wave 
vector QpDW = ^- The CDW order parameters are gen- 
erally incommensurate and have ordering wave vectors 
QcDW — 'z/2,7r ± q/2 (or, equivalently /c2 — ki, 2ki and 

2fc2). 

We will now proceed to write down the bosonized ex- 
pressions of the SC and CDW order parameters. The 
bosonized expressions of the SC order parameters are 



A+ oc e '^'''=+{cos(V7^6'c-)cos(^/7^(/)5+)cos(V7^(/;'s-) 

-I- isin(\/7r6'c-) sin(0r0s+) sin(0r(?is_)} 

(6.20a) 

A_ CX e~'^^<=+ { CO8{y/TT0c-) Sm{y/T:(t)s+) Sm{y/TT(t)s-) 

+ i sin(\/7r6'c-) cos{y/TT(f>s+) cos(\/7r0s-)} 

(6.20b) 

A+ CX 6"'^^"=+ { - COs{^/n(j)c-) COs{y/n(j)s+) COs{y/TT9s-) 

+ isin(\/7r0c-) sin(\/7r(/'s+) sin(\/7r^s-)} 

(6.20c) 

A„ CX 6"'^^"=+ { cos(V7F0c-) sin(y^0s+) sm{y/^ds^) 

— i sin(-\/7r(/!)c-) cos{y/TT(j)s+) cos{y/n6s-)} 

(6.20d) 



Here, and below, in order to simplify the notation we have 
dropped the prefactors of these expressions, including the 
Klein factors, whose effects are taken into account in our 
results. (A discussion of the role of Klein factors in the 
identification of phases in ladders is found in ref . [15] . 

The bosonized form of the CDW order parameters n± 



n+ CX e '\/^'^<:+| — cos(\/7r(/'c-) cos(-v/7r0s+) cos(\/7r'/'s-) 

-|- isin(v'7r0c-) sin(-\/7r0s+) sin{y/TT(f>s^)} 

(6.21a) 

n_ oc e~'^'^'=+ { cos{\/tt (j>c~) siii{^/Tr (l)s+) sin(-\/7r(/)s_) 

— isin(-\/7r?!'c-) cos{^/Tr(l)s+) cos{\/TT<j)s-)} 

(6.21b) 

n+ oc e~'^'^'=+ { — cos(^/7r6'c-) cos(V7r(/>s-i-) cos(-/7r0s-) 
+ isin{y/7r6c-) sin{y/n(l)s+) sin{y/'K6s_)} 

(6.21c) 

n_ oc 6"'^'^'=+ { cos(^/^6'c-) sin(v^0^+) sin(^/^6',_) 

— isin(v'vr6'c-) cos(A/7r0s+) cos(\/7r6's-)} 

(6.21d) 

where we have also dropped the prefactors and their 
dependence on the Klein factors. 



The effective field theory of Eq.(6.12) shows that the 
spin sector s+ couples to the two remaining sectors, the 
charge sector c— and the spin sector s— , only through 
terms that involve the operator cos(-\/47r</)s4_) but not 
the dual field 9s+- A consequence of this feature of 
the effective Hamiltonian is that the Luttinger param- 
eter Ks+ always decreases under the RG flow, as can 
be seen by an examination of Eq.(Alb|, and flows to a 
regime in which Ks+ 0. In this regime the field (ps+ 
is locked and its fluctuations become massive. Hence 
there is a gap in the spin sector, the fleld (f>s+ is pinned, 
and (cos(-\/47r(/)s+)) ^ has a non- vanishing expectation 
value. 

The RG equations given in Appendix [A] reveal that 
for the range of parameters of physical interest all the 
coupling constants (including those in Eq.( |6.22| ) generi- 
cally flow to strong coupling. Hence, we expect that the 
operators cos(\/47r(?!'s-), cos(\/47r0s_), cos(-\/47r0c-), and 
cos(-\/47r0c-) will acquire an expectation value and that 
the fields become locked to the values (/>c- — n^^_y/TT/2, 
Gc- = n0^_^/2, 4>s_ = n^^_0r/2, where n^^_, ne^_, 
ri0^_ , and nff^_ are integers that can each be even or 
odd. Depending of this choice the locked states repre- 
sent different phases. In addition, we recall that oper- 
ators involving dual fields cannot have an expectation 
value simultaneously as this is forbidden by the commu- 
tation relations. This leads us to the conclusion that 
in general we will have different phases depending on 
which fields are locked and to which values. We will 
label the phases by the locked fields: {(f>c-,(t>s-,4's+), 
{(j)c-,Os-,cj)s+), {0c~,(t>s^,<f>s+), and {9c-,9s-,(f)s+) re- 
spectively. Thus, in general we will have a total of eight 
phases characterized by different order parameters. In 
all these phases only the charge sector c+ remains gap- 
less. Additional gapless excitations appear at the contin- 
uous quantum phase transitions between these different 
phases. 

From the structure of the effective field theory we see 
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that the c+ charge sector decouples and remains criti- 
cal for all values of the parameters. It is an effective 
Luttinger liquid with Luttinger parameter Kc+ and ve- 
locity Vc+ ■ This sector has the trivially self-duality of the 
Luttinger models, which guarantees the existence in the 
phase diagram of a dual CDW state for any SC state, 
and vice versa. We will denote this duality symmetry by 



a. 



Uniform SC phases: The bosonized expressions 



of Eq.(6.20a) and Eq.(6.20b) for the two uniform SC or- 



der parameters, A±, imply that these operators may ex- 
hibit quasi long range order provided that the c— sec- 
tor is gapped such that the dual field 9^- is pinned and 
its vertex operator cos^^/ttOc-) has a nonzero expecta- 
tion value. Thus, the uniform SC A_|_ phase (even under 
the exchange of the two legs) occurs whenever the fields 
lock to the classical values {9c-,(j)s~-i4's+) = (0,0,0) or 
{9c-t4's-t4's+) = (""/S, 7r/2, 7r/2). Similarly, the uni- 
form SC A_ phase (odd under the exchange of the two 
legs) occurs whenever the fields lock to the classical val- 
ues (61c-, </>«-,(?!)«+) = (0,7r/2,7r/2) or {9c-As-,(i)s+) = 
(V2,0,0). 

h. PDW phases: The PDW phase A4. occurs for 
{<^,-,9,-,<j)s+) = (0,0,0) and (0c-A-,0s+) = 
(7r/2, 7r/2, 7r/2), while the PDW phase A_ occurs for 
(0c-,0s-,<^s+) = (0,7r/2,7r/2) and (0^-, 0.-, </'s+) = 
(7r/2, 0,0). As it should, the order parameters A-t and 
A±, which describe PDW phases which are even and odd 
under the exchange of the two legs respectively, exhibit 
power law correlations due to the contributions form the 
charge c+ sector. Comparing the bosonized expressions 
for A± and A-t it is clear that uniform SC phases and 
PDW phases are related by the combined dual transfor- 
mation of the two sectors, x 1,?^ . In this system 
PDW phases cannot occur in the absence of Umklapp 
process available at flux $ = tt, and for this reason are 
absent for other values of the flux. 

c. CDW phases: Similarly, the CDW phase 
rij^ has quasi long range order if the field that 
now lock are (0c-j 0s-, 0s+) — (0,0,0) or 
(0c-,0s-,0s+) = {-n 12,71/2,71/2), the phase n_ for 
(0c-,0s-,0s+) = (0,7r/2,7r/2) or (0^-, 0s-, 0s+) 
(7r/2,0,0), the phase n+ for {9 c- , 9 s- , (j) s+) = (0,0,0) 
or {9c-,9s-,4>s+) — (7r/2, 7r/2, 7r/2), and n- for 
(0e-,^s-,0s+) = (0,V2,7r/2) or (^c-, f?s-, 0s+) = 
(V2,0,0). 

The diagram of FigjS] illustrates the symmetry rela- 
tions between various order parameters. 



C. Quantum Phase Transitions 

The effective field theory of the ladder with flux $ = tt 
given in Eq.( 6.12 1 has many effective parameters and cou- 
pling constants. We will not attempt to give a detailed 
description of this theory here. Some important details 
are given in the Appendices. In particular the RG equa- 
tions for the effective field theory are given in Appendix 
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FIG. 5. The relation between various uniform and staggered 
SC and their CDW counterparts present in the phase diagram 
of the model with the flux ^ — tt per plaquette. 



[A| and their solution for general couplings is a complex 
problem. From the simpler case of the standard ladder we 
know that there is always a regime in which the couplings 
flow to strong values which in this case also corresponds 
to a system with a spin gap. The situation is very sim- 
ilar here. Thus while there are regimes in which some 
sectors can remain gapless, there is also a generic regime 
in which only one sector, the charge c+ sector, remains 
gapless while all the other ones are massive. 

Let us now look at the effective field theory under the 
assumption that the s+ sector is massive (and hence 0s+ 
is pinned). We will now examine in detail the dynamics 
of the two remaining sectors, the charge sector c— and 
the spin sector s— . In this regime Eq.(6.12| reduces to 



the simpler system (ignoring for now the decoupled and 
critical charge sector c+) 



H 



Vc- 

2 

V 



+ ^ {K,^{d,9s-f + x;i(a,0,-)2} 

-I- g*i cos(V47r0s-) -I- g*2 cos{V^9s-) 
+ gl cos(V47r0c-) + glz cos(V47r6'c-) 
cos(-\/47r0c-) 



2(7ra)2 

cos(\/47r0c-) " 
2(7ra)2 



33 cos(\/47r6's-) +54 cos(\/47r0s-) 



5„3 cos(V47r6's-) -f- .g„4 cos(V47r0s-) 

(6.22) 



where we absorbed the expectation values of the 
s-\- sector in the effective coupling constants — 
2gQ(cos(%/47r0s+))/(27ra)2, where ga = ffsi, 5s2, ffs, S^s 
respectively. 

Let us consider the subspace defined by gs2 — gs ~ 
5«3 = in the parameter space. From the RG equa- 
tions it can be inferred that once we start with the this 
initial condition, the RG fiow will remain on the same 
hypersurface defined by (7^2 = .93 = Qus = 0. In this 
regime luttinger parameters Kc± and Ks± follow the RG 
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equations below 



VII. CONCLUSIONS 



dK, 



c+ 



dl 
dKs+ 

dl 
dKs- 

dl 
dKc- 

dl 



= 



^^s+ , 2 I 2\ 
-^(5.1+55) 

'^(.9.1 +.94) 



87r2 



(6.23a) 
(6.23b) 
(6.23c) 
(6.23ci) 



The first equation states that the Luttinger parameter of 
the decoupled c+ sector does not renormalize. The sec- 
ond and the third equations state that Ks± renormalize 
to small values, Ks± — > 0. This means both s± sectors 
are opening up a gap such that the vertex functions of 
(j}s± will acquire expectation values. The effective Hamil- 
tonian when the for the regime where (j)s+ is pinned is 



n 



+ ^ cos(%/4^6l,_) + COs(N/4^(/)e_) 



(6.24) 



in which 



where Cg 



9b - §¥^54 and 5^ = §f^gu4 wiieie 

(cos(-\/47r0s+)). This is the same effective theory of 
Eq. ( |4.11 1 in section IV E except that it is written for the 
dual fields in the charge c— sector (instead of the spin sec- 
tor) . It predicts existence of a pair of dual phases which 
between them there is a phase transition in Ising critical 
class. The duality symmetry in the ( 6.24[ ) will be denoted 
by ^2". It relates the state presented by SC operators 
A± to the states with the same parity presented by n±. 
Similarly A-t phases and h± are dual under x Z^~. 
Similar analysis holds in the s— sector. The states with 
the same parity in (A±, n±) and (A-t, n±) are dual under 
Zf- X Z^-. 

On the other hand, if we assume that there is relation 
between some of the couplings (up to restrictions imposed 
by the SU{2) spin invariance) we arrive to a system that 
can be solved by refcrmionization. This is discussed in 
detail in Appendix |B] Depending on the relations be- 
tween the coupling constants the system may be in one 
of the phases we discussed above or be qunatum critical. 
We find two types of quantum criticality. One possibility 
is an ising quantum critical point at which one of the Ma- 
jorana fermions becomes massless. Clearly we have four 
choices for this. On the other hand we also find a case in 
which two Majorana fermions become massless. In this 
case the system has a quantum critical regime which can 
be described as an effective Luttinger model coupled to 
a massive Thirring model. Away from the quantum crit- 
ical regime this system becomes a theory of four coupled 
massive Majorana fermions. 



In this paper we investigated the mechanisms of forma- 
tion of pair-density- wave superconducting order in quasi- 
one-dimensional systems. Although at present time the 
existence and relevance of the PDW state to the inter- 
pretation of experiments in the cuprate superconductors 
can be argued on purely phenomenological grounds, we 
know that this is not a state that is naturally favored 
in a weak coupling BCS theory. The main motivation of 
this work is to investigate the mechanisms of formation of 
PDW order. For this reason it is natural to examine how 
(and if) it it appears on one and quasi-one-dimensional 
systems. 

Here we investigate the occurrence of PDW phases in 
two models of two-leg ladders. In the first model we re- 
examined the properties of the spin-gap phase of a model 
of a two-leg ladder in the regime where the microscopic 
interactions are repulsive and showed that it includes a 
phase with PDW order. Here we showed that within the 
repulsive regime, a PDW state exists provided that one 
of the bands, the bonding band for example, is kept at 
half filling. We showed that in this regime the phase dia- 
gram of the ladder has, in addition to a conventional Lut- 
tinger liquid phase, two superconducting phases: a phase 
with uniform superconducting order (with power law cor- 
relations) and a PDW phase, a superconducting state 
(again with power law correlations) but with wave vector 
QpDW = ^- We also investigated the nature of the quan- 
tum phase transition between these two superconduct- 
ing states and showed that it is in the Ising universality 
class. We discussed in detail the connections that exist 
between this system and the Kondo-Heisenberg chain. In 
particular, much as in the case of the Kondo-Heisenberg 
chain, the PDW order parameter in the two-leg ladder 
is a composite operators of two order parameters of the 
bonding and anti-bonding bands which separately have 
only short range order. Thus this is a highly non-BCS 
realization of PDW order. By extending the analysis to 
the case other commensurate fillings of the bonding band, 
we showed that the state with PDW order arises in con- 
junction with the development of a commensurate CDW 
state. In this sense this result embodies the notion of 
intertwined orders proposed in Ref . ^ . 

We also investigated the existence of PDW phases in 
an extended Hubbard-Heisenberg model on a two leg lad- 
der with flux $ per plaquette. We showed that com- 
mensurate PDW phases appears in this system when the 
flux $ = TT per plaquette. In contrast to the case of 
the conventional ladder, this realization of PDW order 
in the flux $ = 1 ladder can be expressed as a bilin- 
ear of fermion operators. In this sense this realization of 
the PDW state is closer in spirit to the construction of 
FFLO states although in the problem at hand the spin 
rotational symmetry is kept unbroken at all levels. PDW 
order also appears at other values of the flux but only 
when certain commensurability conditions are met, just 
as it is the case in the conventional two-leg ladder. 
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There are still several interesting open questions. 
While the results of this work, and the earlier results 
of Ref. [TT] , show how the pair-density- wave state arises 
together with a spin gap in a system with repulsive in- 
teractions, the ordering wave vector we find is always 
commensurate. However there is no reason of principle 
for the PDW ordering wave vector to be commensurate. 
The root of this phenomenon is the magnetic mechanism 
of the PDW order which is present in both the two-leg 
ladder and in the Kondo-Heisenberg chain. Indeed in 
both cases the ordering wave vectors of the PDW and of 
the spin order (even though it becomes short ranged by 
the development of the spin gap) are the same. On the 
other hand, it is not possible to have incommensurate 
magnetic order (even with power law correlations) in one 
dimension with full SU{2) spin rotational invariance. In- 
deed it is known from work in frustrated one-dimensional 
systems that the incommensurate magnetic state is pre- 
empted in one dimension by a dimerized state with a spin 
gap. Naturally, one way around this problem is to con- 
sider systems with a weak magnetic anisotropy. At any 
rate the construction of a system with incommensurate 
PDW order is an interesting open problem. 



Appendix A: RG equations for the flux ^ — tt model 



The RG equations for the model with flux $ = tt per 
plaquette are 



dK, 



c+ 



dl 

dKs+ 
dl 

dKr- 



8^2 (9sl 



2 I 2 I 2 \ 
?s2 +55 +9u5) 



(Ala) 
(Alb) 



1 

f 2 , 2 I 2\ "c- / 2 I 2 I 2 \ 
-^{93+94+95)- [9u3 + 9u4 + 9u5) 



dl 



dK,^ 



dl 



87r2 



dgs 



dl 
dgs2 
dl 

dl 
dg4 
dl 
dg5 
dl 

dgu3 
dl 

dgu4 

dl 
dgub 

dl 



(Ale) 

(.9s 1 + 54 + 3«4) + ^ (5s2 + .93 + 9I3) 

(Aid) 

(Ale) 



{2-Ks+-K,^)gsi 
1 , 



(2 - Ks 



(2- 
(2- 
(2- 



1 



1 

1 
1 



)9s2 

).93 - 



1 

8V2 

.9455 9u49u5 

~2n 27r 

.93.95 9u39u5 



2tt 



(2 - i^c- 

(2 - i^c- 
(2 - K,^ 



Ks-)94- 
Ks+)95 - 



1 



■)9u3 



Ks-)gu4 



27r 

9s295 

9 si 9 5 
2tt 

9 si 94 _ 5^253 
27r 2tt 

_ 9s29u5 

27r 

_ 9sl9u5 

2tt 



(Alf) 
(Alg) 
(Alh) 
(Ali) 

(Alj) 
(Alk) 



\ 9si9u4 9s29u3 /Ain 
I^s+)9u5 ^ ^ (All) 



2tt 



2tt 



ACKNOWLEDGMENTS 



with the extra constraint (75(0) = 54(0) + 33(0) and 
5u5(0) = 5«4(0) + gu3{0) to guarantee the spin SU{2) 
symmetry. The above set of RG equations, just as the 
Hamiltonian itself, is invariant under all the duality sym- 
metries defined in Section IVII as well as under the ex- 
change of the two Fermi points, 1 O 2. 



We thank Erez Berg and Steven Kivelson for very stim- 
ulating discussions and a previous collaboration which 
motivated this work, and G. Roux for point us out Refs. 
|18' and '39', This work was supported in part by the Na- 
tional Science Foundation, under grants DMR 0758462 
and DMR-1064319 (EF) at the University of lUinois, and 
by the U.S. Department of Energy, Division of Mate- 
rials Sciences under Award No. DE-FG02-07ER46453 
through the Frederick Seitz Materials Research Labora- 
tory of the University of Illinois. 



Appendix B: Refer mionized effective field theory for 
the two-leg ladder with flux $ = tt 



Here we will assume that the s+ sector is gapped and 
the the gapless charge sector c+ is decoupled. The ef- 
fective Hamiltonian for the coupled c— and s— sectors is 
given in Eq.(6.22). Here we will discuss the refermion- 



ized version of this effective field theory for some special 
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combinations of parameters. 



+ 



Cs+ 

2(7ra)2 



gsi cos(\/ 47r03_) + 5^2 cos(V47r6'5-) 



+ 55 cos(v47r0c-) +5«5 cos(v47re'c-) 
cos(\/47r6'c- 



2(7ra)2 

cos(\/47r(/)c_" 
2(7ra)2 



(73 cos(v47r6's-) + 54 cos(v47r0s_) 



5i„3 cos(V47r6's_) +g„4 cos(v47r(/)5_) 



(Bl) 



where Cs+ ~ (cos(v47r(/)s+)). Spin SU{2) invariance 
dictates 33 + 54 = 55 and g„3 + gu4: = g«5- Let us 
assume that — gu3 — gi = gui = 5/2 which im- 
plies that 55 = 5,(5 = g. Moreover let's assume that 
gsi = gs2 = 55 = 5u5 = 5, which implies Ks-{0) = 1. 
For this regime of parameters , and assuming that the 
velocities Vc- = Vg- — v are equal, the Hamiltonian sim- 
plifies to (with a ~ c— , s— ) 



2(7ra)2 



{ cos(V47r0Q) + cos{V4jr9a)} + gl++ 



(B2) 



where the fields I^a' are defined as 



(27ra)2 



cos(v47r0c-) + cr cos (v47r^?c-) 
cos(\/47r0s_) -I- a' cos{\/A^cf)s-) 



(B3) 

We will now assume that also K^- = 1. (Below we will 
relax this assumption.) We can now define two species 
of chiral Majorana fermions 



x2 + ie2 = 



t/4 



7ra 

t/4 



-iv 4:7r{pL^^ 



(B4) 



where a = c— , s— . 

It can be shown that the Majorana mass terms have 
the bosonized form 



mi 



1 

27ra 
1 

2'Ka 



cos(v 47r0ct) -I- cos( V 47r0Q) 

COs{V^(t)a) - C0s{V^9a) (B5) 



Using the equations above, the total Hamiltonian for this 
sector reads as (after setting the velocity v = 1) 



n = - 



E 

a — c—.s — 



E 
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(B6) 



where the Majorana mass is M = gCs+fiTra). Here, 
(^fl '^L ) ^^'^ (^iT'^L ) massless Majorana 

fields. In this case the system is at a quantum critical 
point. 

We now note that if Kc- 7^ 1 is allowed, the 
refermionized theory now has a Luttinger-Thirring four 
fermion (current current) coupling term for the fermions 
in the c— sector of the form gR}^_Rc-L\_Lc- = 
^4.9 Xfl~'^fl~Xir^L~) where g measures the departure 
from Kc- = 1. This term, which in the conventional 
Luttinger-Thirring model is marginal, in this case mixes 
the massless sector with the massive sector. However if 
we were to integrate out the massive sector it will induce 
a marginal operator in the remaining massless fermions. 
We will see below that the same marginal operator is 
present automatically if we relax some of the relations 
between the coupling constants. For this reason we will 
ignore these terms for the time being. 

On the other hand, if some of the relations between 
the coupling constants are lifted (but keeping track of 
the constraints due to the SU (2) spin symmetry) all four 
Majorana fields become separately massive and the sys- 
tem is in one of the phases described in section [VI B[ In 
this language we can picture the system becoming quan- 
tum critical by turning one or two Majorana fermions 
massless. The case with one Majorana fermion becom- 
ing massless is the Ising quantum criticality that we have 
already discussed. 

Let us now focus on the case in which one pair of 
Majorana fields remains massless. In this case we can 
build Dirac Fermions out of the Majorana fermions. This 
transformation will mix the charge and spin fields into a 
new Dirac field. The right- and the left-moving compo- 
nents are defined as 



(B7) 



In terms of the new variables, the Hamiltonian reads as 
^ = E + MMi + f (^17'^V'i)' (B8) 

i=l,2 



where = cTx, 1^ = — io'yj 7^ — 1^1^ = fz and 



ipj = {RI,lI). This Hamiltonian consists of a free 
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fermion sector described by free Dirac field V'2 while the 
dynamics of the ipi is described by massive Thirring 
model. Therefore in this regime of parameters the total 
symmetry associated with the s— and c— sectors is 
U{1) X U{1) which includes a global C/(l) symmetry and 
the chiral symmetry of the free sector generated by 75. 

For the regime of parameters ^si = 5ti5 = 5i and gs2 — 
.95 = .92 while gsi 7^ 5s2, the ^2 Fermion becomes massive. 
Thus, the effective Hamiltonian of Eq.(B8) now has two 
mass terms of the form 



(B9) 



with masses given by and M = -^q+ and m = -^q^ 
where g± = {gsi ± 3s2)/2. In this language we have four 
phases depending on the relative signs of the masses M 
and m of the two Dirac fermions. 

Under the these assumptions, except for the SU{2) in- 
variance, the {53,54,(7113,51,4} are arbitrary. We re-write 
the four-fermion interactions as 



— fii^++ + fi2^A — + ./21I- 



/222:- 



(BIO) 



where new coupling constants /o-a' is related to the 5's 
as 





(Bll) 



Using Eq. ( B5 ) , one can write lo-o-' in terms of the Ma- 
jorana fermions. as we saw earlier, is just the 

current-current interaction of the form (V'iT'^V'i)^- By 

similar argument, I is the same type of interaction 

but for the -02 fields. Together they add up to the fol- 
lowing 

fia+++f22^— = /ll(^l7^V'l)'+/22(^27^^2)' (B12) 

where ip = tp^jo — {L\R}). Each operator Xj^+ and 
T are invariant under the t/(l) x U(l) symmetry com- 
posed of a global U{1) symmetry associated with charge 
conservation and the continuous chiral symmetry of each 



Dirac fermion. The mass terms break the chiral symme- 
try down to a discrete Z2 symmetry. 

However, the off-diagonal terms I+- and T ^ involve 

both X and f fields. In terms of Majorana fermions they 
read 



(B13) 



As it turns out these terms violate the conservation of 
fermion number of the Dirac fermions and, for this rea- 
son, are more naturally expressed in terms of Majorana 
fields. 

Let us look at the regime of parameters in which the 
Dirac fermion number violating couplings are absent and 
set /12 = /21 — 0. This happens when 53 = 5^4 and 54 = 
5„3. According to the SU{2) spin rotation invariance 
condition, in this regime 55 = 5^5, and therefore the 
Dirac fermion ipi remains massless. Assuming gsi = gs2 
the Hamiltonian will read as 



Gi 



G2 



(V527^V'2)' 



(B14) 



where Gi = /n =33-1-34 and G2 = /22 = 53 - 54- 
The resulting Hamiltonian splits into the Hamiltonian of 
massless Thirring model for il)2 and a massive Thirring 
model for the Tpi with mass M = ^^(53 -I- 54). There- 
fore what has changed with respect to the case when 
53 = 54 = 5«3 = 5m4 is that the dynamics of V'2 
field is now described by the mass-less Thirring model 
which shares the same U{\) x U{\) symmetry with the 
non-interacting case G2 — 0. Thus, under the (less - 
restrictive) conditions 53 — 5^4 and 54 — g^s, the sym- 
metry is still U{1) X [/(I). 

Therefore in this case the system decouples into a 
massless Thirring model and a massive Thirring model 
each with a separate conserved charge current. The mas- 
sive Thirring model is an integrable system which by 
bosonization can be mapped onto the sine Gordon field 
theoryjSsl in the regime in which the sine Gordon term 
is relevant. Hence this sector has a spectrum of mas- 
sive solitons. On the other hand, the massless Thirring 
model, which is equivalent to a spinless Luttinger model, 
is a quantum critical system with an exactly marginal op- 
erator, parametrized by the coupling constant G2- Hence 
in this case instead of Ising quantum criticality we get 
Luttinger quantum criticality. 
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